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Abstract 

' We consider Brownian motion in a bounded domain $7 on a two- 

£SJ ■ dimensional Riemannian manifold (S, <?). We assume that the boundary 

d£l is smooth and reflects the trajectories, except for a small absorbing 
arc d£l a C dCl. As <9il a is shrunk to zero the expected time to absorption 
in dQ a becomes infinite. The narrow escape problem consists in con- 
O structing an asymptotic expansion of the expected lifetime, denoted Et, 

as e = |<9r2 a | 9 /|<9r2| 9 — > 0. We derive a leading order asymptotic approx- 

1 1 

log — h 0(1) . The order 1 term can be evaluated 

for simply connected domains on a sphere by projecting stereographi- 

cally on the complex plane and mapping conformally on a circular disk. 

^ I It can also be evaluated for domains that can be mapped conformally 

$— i ' onto an annulus. This term is needed in real life applications, such as 

1 

trafficking of receptors on neuronal spines, because log — is not neces- 

£ 



d • imation Et = — - 

Dtt 



sarily large, even when e is small. If the absorbing window is located 



at a corner of angle a, then Et — 



log i + o(i) 



if near a cusp, 



Da 

then Et grows algebraically, rather than logarithmically. Thus, in the 
domain bounded between two tangent circles, the expected lifetime is 
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1 Introduction 



In many applications it is necessary to find the mean first passage time (MFPT) 
of a Brownian particle to a small absorbing window in the otherwise reflecting 
boundary of a given bounded domain. This is the case, for example, in the 
permeation of ions through protein channels of cell membranes JI], and in the 
trafficking of AMPA receptors on nerve cell membranes [2] , [H] • While the first 
example is three dimensional the second is two dimensional, which leads to 
very different results. In this paper we consider the two dimensional case. 

In the first two parts of this series of papers, we considered the narrow 
escape problem in three dimensions J3J and in the planar circular disk j3] . The 
leading order asymptotic behavior of the MFPT is different in the three and 
two dimensional cases; it is proportional to the relative size of the reflecting 
and absorbing boundaries in three dimensions, but in two dimensions it is 
proportional to the logarithm of this quotient. The difference in the orders of 
magnitude is the result of the different singularities of Neumann's function for 
Laplace's equation in the two cases. 

While the second term in the asymptotic expansion of the MFPT in three 
dimensions is much smaller then the first one, it is not necessarily so in two 
dimensions, because of the slow growth of the logarithmic function. It is 
necessary, therefore, to find the second term in the expansion in the two- 
dimensional case. This term was found for the case of a planar circular disk in 
jS], and can therefore be found for all simply connected domains in the plane 
that can be mapped conformally onto the disk. Similarly, it can be found for 
simply connected domains on two-dimensional Riemannian manifolds that can 
be mapped conformally on the planar disk. For example, the sphere with a 
circular cap cut off can be projected stereographically onto the disk, and so 
the second term for the narrow escape problem for such domains can be found. 

The specific mathematical problem can be formulated as follows. A Brow- 
nian particle diffuses freely in a bounded domain flona two-dimensional Rie- 
mannian manifold (£, g). The boundary dfl is reflecting, except for a small 
absorbing arc dfl a . The ratio between the arclength of the absorbing boundary 
and the arclength of the entire boundary is a small parameter, 

\dn\ 9 * 

The MFPT to dfl a , denoted Et, becomes infinite as e — > 0. 

In this paper we calculate the first term in the asymptotic expansion of 
Er for a general smooth bounded domain on a general two-dimensional Rie- 
mannian manifold. We find the second term for an annulus of two concentric 
circles, with a small hole located on its inner boundary. This result is general- 
ized in a straightforward manner to domains that are conformally equivalent 
to the annulus. 
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The calculation of the second term involves the solution of the mixed 
Dirichlet-Neumann problem for harmonic functions in Q. While in the three 
dimensional case this is a classical problem in mechanics, diffusion, elastic- 
ity theory, hydrodynamics, and electrostatics the two dimensional 
problem did not draw as much attention in the literature. 

First, we consider the problem of narrow escape on two dimensional man- 
ifolds, and derive the leading order asymptotic approximation 



Et 



Dtt 



io g i + o(i; 



for £<1. (1.1} 



This generalizes the result of [2 J from general smooth planar domains to general 
domains on general smooth two-dimensional Riemannian manifolds. 

The second term in the asymptotic expansion is found for the 2-sphere 
x 2 + y 2 + z 2 = R 2 . The calculation is made possible by the stereographic 
projection that maps the Riemann sphere onto a circular disk, a problem that 
was solved in [Hj. The boundary in this case is a spherical cap of central angle 
5 at the north pole, where e is the ratio between the absorbing arc and the 
entire boundary circle. We find that the MFPT, averaged with respect to an 
initial uniform distribution, is given by 



Et 
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2nD 



log - + 2 log - + 3 log 2 - I + 0(e, 6 2 log 6, 5 2 log e) 
o e I 



(1.2) 



where |Q| fl = 4ttR is the surface area of the sphere. Note that there are two 
small parameters that control the behavior of the MFPT in this problem. The 
small e contributes as equation (jl.lj) predicts, whereas the small 5 parameter 
contributes half as much. 

The second case that we consider is that of narrow escape from an annulus, 
whose boundary is reflecting, except for a small absorbing arc on the inner 
circle. Specifically, the annulus is the domain Ri < r < R 2 , with all reflecting 
boundaries except for a small absorbing window located at the inner circle (see 
Fig. P). The inversion w = 1/z transforms this case into that of the absorbing 

R 

boundary on the outer circle. Setting (3 = — < 1, the MFPT, averaged with 

R2 

respect to a uniform initial distribution, can be written as 
Et = 



CP 2 p 2n 



log- + log2 + 2/? 2 

e 



1 

+ ^ 



R\ 



2 1 - P 2 /3 



log---R 2 2 + 0(e,[3 4 )R 2 2 



Also in this case we find two small parameters, the e contribution belongs to a 
singular perturbation problem with a boundary layer solution and an almost 
constant outer solution with singular fluxes near the edges of the window, 
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whereas the 5 contribution is just the singularity of Green's function at the 
origin-a problem with a regular flux. This result is generalized to a sphere with 
two antipodal circular caps removed. We find that for [3 <C 1 the maximum exit 
time is attained near the south pole, as expected. This result can be generalized 
to manifolds that can be mapped conformally onto the said domain. 

The asymptotic expansion of the MFPT to a non-smooth part of the bound- 
ary is different. We consider two types of singular boundary points: corners 
and cusps. If the absorbing arc is located at a corner of angle a, the MFPT is 



Et 



\Q 
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log i + 0(1) 



For example, the MFPT from a rectangle with sides a and b to an absorbing 
window of size e at the corner (a = 7r/2, see Figure |2J), is 

a 2 7i b 



2fi 

Et = -!— ^ 

7T 



log- + log- + ^ + 2/3 2 + 0^,/3 4 ) 
e tt b a \a J 



where |0| = ab and (3 = e nb ^ a . The calculation of the second order term turns 

out to be similar to that in the annulus case. The pre-logarithmic factor is 

Da 

the result of the different singularity of the Neumann function at the corner. It 
can be obtained by either the method of images, or by the conformal mapping 
z i — > z n / a that flattens the corner. In the vicinity of a cusp a — > 0, therefore 
the asymptotic expansion (jl.3j) is invalid. We find that near a cusp the MFPT 

grows algebraically fast as — r , where A is the order of the cusp. Note that 
the MFPT grows faster to infinity as the boundary is more singular. The 
change of behavior from a logarithmic growth to an algebraic one expresses 
the fact that entering a cusp is a rare Brownian event. For example, the 
MFPT from the domain bounded between two tangent circles to a small arc 

at the common point (see Figure EJ) is Et = j — — r — ( — h 0(1) V where 

[a 1 — l)D \e J 

d < 1 is the ratio of the radii. This result is obtained by mapping the cusped 
domain conformally onto the upper half plane. The singularity of the Neumann 
function is transformed as well. The leading order term of the asymptotics can 
be found for any domain that can be mapped conformally to the upper half 
plane. 

In three dimensions the class of isolated singularities of the boundary is 
much richer than in the plane. The results of [3] cannot be generalized in a 
straightforward way to windows located near a singular point or arc of the 
boundary. We postpone the investigation of the MFPT to windows at isolated 
singular points in three dimensions to a future paper. 

As a possible application of the present results, we mention the calculation 
of the diffusion coefficient from the statistics of the lifetime of a receptor in a 
corral on the surface of a neuronal spine 



4 



2 Asymptotic approximation to the MFPT on 
a Riemannian manifold 



We denote by x(t) the trajectory of a Brownian motion in a bounded domain Q 
on a two-dimensional Riemannian manifold (E, g). For a domain Q C E with a 
smooth boundary dfl (at least C 1 ), we denote by the Riemannian surface 
area of Q and by \dfl\ g the arclength of its boundary, computed with respect 
to the metric g. The boundary dQ is partitioned into an absorbing arc dfl a 
and the remaining part dQ — dQ a is reflecting for the Brownian trajectories. 
We assume that the absorbing part is small, that is, 



\dn a \ 
ion 



< i, 
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however, E and Q are independent of e; only the partition of the boundary dfl 
into absorbing and reflecting parts varies with e. 

The first passage time r of the Brownian motion from Q to dQ a has a finite 
mean and we define 

u{x) = E[t\x(0) = x\. 

The function u(x) satisfies the mixed Neumann-Dirichlet boundary value prob- 
lem (see for example P, 

DA g u(x) = -1 for x eQ (2.1) 



du(x) 



for x G dVt — dtt a (2.2) 



dn 

u(x) = for xe dVl a , (2.3) 

where A g is the Laplace-Beltrami operator on E and D is the diffusion coeffi- 
cient. Obviously, u(x) — > oo as e — > 0, except for s in a boundary layer near 

2.1 Expression of the MFPT using the Neumann func- 
tion 



We consider the Neumann function defined on E by 

1 

dN(x, y 



A g N(x,y) = -5(x-y) + -—, for x.yeVl (2.4) 



= 0, for x E <9f2, y G f2. 
an 
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The Neumann function N(x, y) is defined up to an additive constant and is 
symmetric [S]. The Neumann function exists for the domain Q, because the 
compatibility condition is satisfied (i.e., both sides of eq. (|2.4j) integrate to 
over Q due to the boundary condition). The Neumann function N(x,y) is 
constructed by using a parametrix H (x, y) , 

H(x,y) = - h{di *> y)) logd(x,y), (2.5) 

where d(x, y) is the Riemannian distance between x and y and h(-) is a regular 
function with compact support, equal to 1 in a neighborhood of y. As a 
consequence of the construction N(x, y) — H(x, y) is a regular function on Q. 

To derive an integral representation of the solution u, we multiply eq. ()2.1|) 
by N(x,y), eq. ([2.4|) by u(x), integrate with respect to x over Q, and use 
Green's formula to obtain the identity 

N(x(S),$) 9u{ * {S)) dS g = --L [ u (x)dV 9 + u($) 



an ® n \®-\g Jn 



[ N(x,£)dV g . (2.6) 
Jn 



The integral 



C £ = — u(x)dV g (2.7) 
Jn 

is an additive constant and the flux on the reflecting boundary vanishes, so we 
rewrite eq. (|2.6|) as 

u(£) = C e + f N(x,$)dV g + f N(x(S),Z) du{ * { n S)) dS g , (2.8) 

where S is the coordinate of a point on dQ a , and dS g is arclength element on 
dfl a associated with the metric g. Setting 

du(x(S)) 
f[b} ~ dn ' 
and choosing ^ G dfl a in eq. ()2.8|) . we obtain 

= C £ + 

Jn Jan a 

The first integral in eq. (j2.8J) is a constant (independent of e), because due to 
the symmetry of N(x, y) eq. ()2.4j) gives the boundary value problem 



/ N{x,£)dV g + [ N(x(S),£)f(S)dS g . (2.9) 
Jn Jan n 



d 



dn{£) 



A f I N{x,i)dV g = for £eft 
q Jn 

N(x,i)dV g = for £G<9Q, 
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whose solution is any constant. Changing the definition of the constant C £ , 
equation (J2.8|) can be written as, 

«(€)=/ N(x(S),£)f(S)dS g + C e , (2.10) 

Jdfla 

and both f(S) and C £ are determined by the absorbing condition ()2.3j) 

0= / N(x(S),Z)f{S)dS g + C 6 for £edn a . (2.11) 
Jdn a 

Equation (j2.11j) has been considered in j2] for a domain E C R 2 as an integral 
equation for f(S) and C £ . 

Actually, the boundary coordinate S can be chosen as arclength on dQ a , 
denoted s. Under the regularity assumptions of the boundary, the normal 
derivative f(s) is a regular function, but develops a singularity as £(s) ap- 
proaches the corner boundary of dQ a in dQ JTU] ■ Both can be determined from 
the representation (J2.1U)) . if all functions in eq. (|2.11|) and the boundary are an- 
alytic. In that case the solution has a series expansion in powers of arclength 
on Q a . The method to compute C £ follows the same step as in [2|. 

2.2 Leading order asymptotics 

Under our assumptions, u(£) — > oo as e — > for any fixed £ G Q, so that 
eq. ([2.7|) implies that C £ —>■ oo as well. It follows from eq. ()2.11j) that the 
integral in (|2.11|) decreases to — oo. 

An origin G dQ a is fixed and the boundary dfl is parameterized by 
(x(s),y(s)). We rescale s so that dQ = \(x(s),y(s)) : — | < s < || and 
( x ( — I) j 2/ ( — I)) = ( x (I) '2/ (I))- We assume that the functions x(s) and 
y(s) are real analytic in the interval 2|s| < 1 and that the absorbing part of 
the boundary dfl a is the arc 

dn a = {(x(s),y(s)) : \s\ < e} . 

The Neumann function can be written as 

N(x,t) = - 1 }-\ogd(x,S)+v N (x,Z), for xeBs(Z), (2.12) 

where Bg(£) is a geodesic ball of radius 5 centered at £ and v n( x ', £) is a regular 
function. We consider a normal geodesic coordinate system (x, y) at the origin, 
such that one of the coordinates coincides with the tangent coordinate to dQ a . 
We choose unit vectors e\, e% as an orthogonal basis in the tangent plane at 
so that for any vector field X = x\e.\ + X2e2, the metric tensor g can be 
written as 

gij = 5ij + e 2 cifjXkXi + o{e 2 ), (2.13) 

kl 
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where \xk\ < 1, because e is small. It follows that for x,y inside the geodesic 
ball or radius e, centered at the origin, d(x, y) = cLe(x, y) + 0(e 2 ), where cLe is 
the Euclidean metric. We can now use the computation given in the Euclidean 
case in j2]. To estimate the solution of equation 1)2.11)1 . we recall that when 
both x and £ are on the boundary, vn(x,£) becomes singular (see [HI p. 247, 
eq.(7.46)]) and the singular part gains a factor of 2, due to the singularity of 
the "image charge". Denoting by the new regular part, equation 1)2.11)1 
becomes 



|s'|<£ L 



v N {x(s );£(s)) 



7T 



f(s')S(ds') = C £ , (2.14) 



where S(ds') is the induced measure element on the boundary, and x = 
(x(s),y(s)), £ = (£(s),r)(s)). Now, we expand the integral in eq. ()2.14jl . as 
in0, 



\ogd(x(s),Z(s')) = log ^( x (s')-as)) 2 + (y(s')- V (s)) 2 (1 + 0(e 2 )) 
and 

oo oo 

S(ds)f(s) = J2fj sjds , v N {x{s')-£{s))S{ds') = J2 v i( s ') sJ ds ' ( 2 - 15 ) 

3=0 j=0 

for \s\ < e, where Vj(s') are known coefficients and fj are unknown coefficients, 
to be determined from eq. 1)2. 14)1 . To expand the logarithmic term in the last 
integral in eq. (|2.14p . we recall that x(s'),y(s'),£(s), and r](s) are analytic func- 
tions of their arguments in the intervals \s\ < e and \s'\ < e, respectively. 
Therefore 

J\s') n log d(x(s),d(s'))ds'= (2.16) 

(sT log V(x(s')-as)) 2 + (y(s')-v(s)) 2 (1 + 0(e 2 )) ds' = 
J\s') n log {\s'-s\ (1 + O ({s' - s) 2 ))} (l + 0(e 2 )) ds'. 

We keep in Taylor's expansion of log {\s' — s\ (1 + O ((s' — s) 2 ))} only the lead- 
ing term, because higher order terms contribute positive powers of e to the 
series 

£ M« - *T ^ = 45 (log. - 1) + 2 J2 (2/ ' n/ (2.17) 
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For even n > 0, we have 



's') n \ok(s - s') 2 ds' 



4 



\oge 



n+l" a (n + 1) 2 
whereas for odd n, we have 



,2i 



-n-2j+l 



- j(n-2j + l)' 



^Tlog^-*') 2 ^' = - 



x s 2j+l £ n-2j 



4^2j + l 



^ 2 j + ln-2j 
Using the above expansion, we rewrite eq. (|2.14j) as 



:2.i8) 



(2.19) 







£ |_Llog [\s' - s\ 2 (1 + O {(s' - S ) 2 )) (1 + 0(, 2 ))] + 



x 



i=o 

and expand in powers of s. At the leading order, we obtain 



e (logs - 1) /o + J] (^^y 10 ^ 



r 2p+l 



(2p + l) 



2 ) 



Uo(s') ds' + O e . 



(2.20) 



Equation (12~2UD and 



1 /•£ 2p+l 



determine the leading order term in the expansion of C £ . Indeed, integrating 
eq. (|2.1j) over the domain Q, we see that the compatibility condition gives 



f(s)S(ds) = 



ai 



(2-21) 



and using the fact that J v (s')S(ds') = 0(e), we find that the leading order 
expansion of C £ in eq. (|2.2(Jj) is 



101 



7T 



log - + 0(1) 



for £<1. 



(2.22) 



If the diffusion coefficient is D, eq. (|2.1U|) gives the MFPT from a point x G Q, 
outside the boundary layer, as 



E[t I x] = u(x) 



nD 



log - + 0(1) 



for £ < 1. 



(2.23) 
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3 The annulus problem 



We consider a Brownian particle that is confined in the annulus R\ < r < R2. 
The particle can exit the annulus through a narrow opening of the inner circle 
(see FigflJ. The MFPT v(x) satisfies 



Av 



— 1, for R t <r < R 2 



(3.1) 



dv 
dr 

dv 
dr 



0, for r = R 2 , 

0, for r = R 1: \9 — ir\ > e, 



v = 

?2 ^2 



0, for r — R\, 



7r < e. 



R 2 -r 1 

The function w = — ^ is a solution of the Dirichlet problem for eq. (j3.1j) 

in the exterior domain of the inner circle r > R\. More specifically, it satisfies 
the boundary value problem 



Aw 



-1, for R 1 <r < R 2 



dw 
dr 

dw 
dr 



-Ri, for r = Ri, 



--i? 2 , for r = R 2 , 



w 



0, for r — R\. 



The function u = v — w satisfies 



Au 



0, for Ri < r < R 2 , 



(3.2) 



du 
dr 

du 
dr 



-R 2 , for r = R 2 



for r = R u 



7T > £, 



U 



0, for r = R t , \9 — tt\ < e. 



Separation of variables produces the solution 



n=l 



cos n6 + a log 



(3.3) 
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where a n , b n and a are to be determined by the boundary conditions. Differ- 
entiating with respect to r yields 



r-| OO 

OU yr-\ 

Or ' 

71=1 



Ro \ R 



n-l 



b n R 2 I R 



n-l 



a 

cos no H — . 

r 



Setting r = i? 2 gives 



n (a n — b n ) cos n# + a 



n=l 



(3.5) 



(3.6) 



therefore, a n = b n and a = ~^R-2-> an d we have 



n=l 



r 



cos nQ H — i? 2 l°g ( tt 
° " R\ 



(3.7) 



The boundary conditions at r = become the dual series equations 

R, ~ 

R2, 



2 



A \fR2Y (- 



cos nO 



n=l 



na r 



R2 

Ri 



n+l 



Ri 

R2 



n-l 



cosn$ 



0, for \$ — 7r| < e 
2R X 



(Rl-Rj), for \e-ir\>e. 



Setting 



Ri 
R2 



R2 
Ri 



n+l 



Ri 

R2 



n-l 



a„, for n > 1, 



and c = ao converts the dual series equations to 



2 ^ 1 



n=l 



cos n8 



0, for 7r — e < 9 < n, 



1 



(3- 



(3.9) 



^nc n cosn9 = -{R 2 2 -R\) } for < 6 < tt - e, (3.10) 

71=1 

2/3 2n 

where if n = — — for n > 1, and ffn = 0, with /3 = — - < 1. Note that 

l + (3 2n ~ R2 

H n = 0(/3 2n ) which tends to zero exponentially fast (much faster than the n _1 
decay required for the Collins method OHH], see also 111)- 

The case H n = was solved in Q. We now try to find the correction of 
that result due to the non vanishing H n . As in jS] the equation 



00 

CO \ ^ Cn 

~2 



n=l 



cos n6 



cos 



hi(t) dt 



2 Je V cos 6 — cos t 
11 



for for < < 71 — e 



defines the function h\{6) uniquely for < 6 < it — e, the coefficients are given 
by 

1 + H r^—t 

c n = — /=-^ / hit) [P n (cost) + P^cost)] dt, (3.11) 
V 2 Jo 

and 

rTT— e 

c = V2 h x {t)dt. (3.12) 
Integrating equation ()3.1()|) gives 

oo 

^c„sinn0= - (Pj-P?)0, for < < tt - e. (3.13) 

n=l 

Substituting eq. (|3.13|) in equation (J3TTTJ) , changing the order of summation 
and integration, while using [TTJ eq. (2.6.31)], 

-= V[P„(cost) + P n _i(cost) sinnO = 2 ' (3.14) 

V2 , v cost — cos & 

v n=l v 

we obtain for < 6 < tt — e, 

f , k f ] dt+ r £ K,(6 1 t)h 1 (t)dt= {Rl ~ R l )e , (3.15) 
Jo v cost — COS V Jo 9™o 



2 cos ■ 



where the kernel K» is 



/3 

oo 



Kp(6,t) = ^-£-^PT„(P n (cost) + P„_i (cost)) sin nfl 

\/2 COS - n=l 

2 

= -2^(1 + cost) sin 6 -l3 2 + 0{l3 4 ). (3.16) 

The infinite sum in eq. ()3.16|) is approximated by its first term, while using the 
first two Legendre polynomials Pq(x) = 1, P\{x) = x. Using Abel's inversion 
formula applied to equation (|3.15|) . we find that 

. u 

h x {t)- K p {t,s)h 1 {s)ds = -± !-/ 2 = du, (3.17) 

Jo 71 dt J v cos u — cost 

where the kernel Kp is 

/. x 1 d /"* if /3 (M,s)sinM 

= "77 / / (3-18) 

^ dt J \/ cos u — cos t 

u 



TT dt Jo ^/ cos u — COS t 
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The substitution 



gives 



cos^-cost (3i9) 



u 



therefore, 



t sm — sin u + 

. 2 = rfn = ^sin 2 -, (3.20) 
V cos u — cos t V2 2 

K p (t,s) = 2/? 2 cos 2 |sint + 0(/? 4 ). (3.21) 

Equation (|3.17|) is a Fredholm integral equation of the second kind for hi, 
of the form 

(I - k(i)h = z, (3.22) 

where 

u 

R 2 -R 2 d f f usm o J 

z(t) = -2 1— / 2 = du. 

7T at Jo v cos it — cos i 

Therefore, we /i can be expanded as 

h = z + i^z + + . . . , (3.23) 

which converges in L 2 . Since Co = V2(h, 1) (eq. (|3.12Jl ). we find an asymptotic 
expansion of the form 



Co 



V2 



(z,l) + (K p z,l) + ... . (3.24) 



The leading order term of this expansion was calculated in j^J. We now esti- 
mate the error term {Kpz, 1), which is also the 0(f3 2 ) correction. Integrating 
by parts and changing the order of integration yields 

u 

Knz{t) = 2(3 2 — -suit / cos 2 -cfe — / 1 =du 

7r J 2 ds Jq vcosm — coss 

d2 d2 /-tt r s u sin — du 



o -Ko — -^1 • / • / 

p sm r / sm s ds I 

TT Jo Jo 



V cos u — COS s 

x/2/? 2 (i? 2 -i? 2 )sint. (3.25) 



Therefore, 



(K p z,l) = V2f3 2 (Rj- Rl) smtdt = 2V2(3 2 (R 2 2 - R 2 ). (3.26) 

Jo 
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We conclude that 



Co = (Rl - R\] 



21og^ + 21og2 + 4/3 2 + 0(e,/^ 



(3.27) 



The MFPT averaged with respect to a uniform initial distribution is 

Co 



Et 



i Ri 



2 + 2 R\ - R{ 



, #2 1 
2 l0 S — ~ 7 JL 2 



(i?2 - R\) 



log 



i?i 4 
log2 + 2/5 2 



(3.28) 



1 



+ 2T^ lo 4~i^ + ° (£ ' m - 



Note that there are two different logarithmic contributions to the MFPT. The 
"narrow escape" small parameter e contributes 



log-. 



7T 



as expected from the general theory (equation ([1.1)) ) 
(3 contributes 

log-. 



2tt 



(3.29) 

whereas the parameter 



(3.30) 



These asymptotics differ by a factor 2, because they account for different sin- 
gular behaviors. The asymptotic expansion (j3.29j) comes out from a singular 
perturbation problem with singular flux near the edges, boundary layer and 
an outer solution, whereas the asymptotics ()3.30|) is an immediate result of the 
singularity of the Neumann function, with a regular flux. 

The exit problem in an annulus with the absorbing window located at 
the outer circle is solved by applying the complex inversion mapping z — > 1/z, 
which maps the annulus into itself and replaces the roles of the inner and outer 
circles. In such case, in the limit /? — 0, the circular disk problem is recovered, 
where it was shown that the maximum exit time is attained at the antipode 
point on the outer circle. Since the inversion mapping exchanges the inner 
and outer circles, we conclude that for the original annulus problem, where 
the absorbing boundary is located at the inner circle, the maximal exit time 
is attained at the inner circle, at the antipode of the center of the hole. This 
point is close by to the hole itself, a result which is somewhat counterintuitive. 
Note that (J3.28)) is valid, with the obvious modifications, for any domain that 
is conformally equivalent to the annulus. 



4 Domains with corners 

Consider a Brownian motion in a rectangle Q = (0, a) x (0, b) of area ab. The 
boundary is reflecting except the small absorbing segment dQ a = [a—e, a] x {b} 
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(see Fig. [2J). The MFPT v(x,y) satisfies the boundary value problem 

Av = — 1, (x, y) G Q, 
v = 0, G <9fi a , 

^ = o, (x, y ) edn-dn a . (4.i) 

The function f = satisfies 

J 2 

Af = -1, (x,y)eQ, 

f = 0, 

^ = 0, (x,y)G{0}x[0,6]U{a}x[0,6]U[0,a]x{0}, 

^ = -6, (z,2/)e[0,a-e]x{&}, (4.2) 

therefore, the function u = v — / satisfies 

Am = 0, (x,y) G f2, 

u = 0, (x,y)edtt a , 
du 

— = 0, (x,y)G{0}x[0,6]U{a}x[0,fe]U[0,a]x{0}, 
(9m 

— = b, (x,y) G [0, a — e] x {b}. (4.3) 
A solution for u in the form of separation of variables is 

oo 

u(x,y) = h > a n cosn cos , (4.4) 

2 / — ' a a 

n=l 

where the coefficients a n are to be determined by the boundary conditions at 
y = b 

oo , 

a o sr^ 7m " 7Tnx \ 

u(x,b) = h > a n cosh cos = 0, x G (a — e,a), 

2 a a 

71=1 

du, , N TTV^ • ! ™ & 7TWX . , 

— (x, o) = — > na n smh cos = o, j;G(0,a-£). (4.5) 

av a ^— f a a 

y n=l 

Setting c n = a n smh , we have 

a 



f 'ii , <",! 

n=l 



2 ^ 1 + H. 



cos n9 = 0, 7r — 5<9<n, 



y^nc n cosn9 = — , < 6 < n - 5, (4.6) 

n=l 
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where 5 



7T5 



and H„ 



Tib 



tanh 



Tinb 



l,n > 1. Note that H n = O (f3 



2n\ 



for (3 = exp < > < 1. The rectangle problem and annulus problem (eq. 

{ a ) 

(13.10)) ) are almost mathematically equivalent, and equation ()3.27j) gives the 
value of Co 



c 



2a6 

7T 
7T 



21og- + 21og2 + 4/3 2 + 0(<5,/? 4 
o 

log - + log - + 2/3 2 + (-,/?' 
en \a 



(4.7) 



The error term due to 0{(3 i ) is generally small. For example, in a square a = b 
and (3 = e~ n so that /3 4 ~ 3 x 10 -6 . The MFPT averaged with respect to a 
uniform initial distribution is 



c 6 2 2a6 

£/T = = 

2 3 7T 



, a i 2 n b ^ n2 (S n . 
log - + log - + -- + 2(3 2 + O - /? 4 
£ 7T o a \a 



The leading order term of the MFPT is 

2|Q 



log 



(4.8) 



(4.9) 



which is twice as large than ()3.29j) . The general result f j 1 . 1 1) was proved for 
a domain with smooth boundary (at least C 1 ). However, in the rectangle 
example, the small hole is located at the corner. The additional factor 2 is the 
result of the different singularity of the Neumann function at the corner, which 
is 4 times larger than that of the Green function. At the corner there are 3 
image charges — the number of images that one sees when standing near two 
perpendicular mirror plates. In general, for a small hole located at a corner of 
an opening angle a (see Fig. EJ), the MFPT is to leading order 



Da 



(4.10) 



7T 

This result is a consequence of the method of images for integer values of — . 

a 

For non-integer — we use the complex mapping z \— > z*l a that flattens the 
a 

corner. The upper half plane Neumann function — log z is mapped to — log z 

___ 7r a 

and the analysis of Section El gives ()4.10j) . 

To see that the area factor \Q\ remains unchanged under the con- 
formal mapping / : (x,y) i— > (u(x,y),v(x,y)), we note that this factor is 
a consequence of the compatibility condition, that relates the area 
to the integral 

A^ y) wdxdy = -- 



J 

Jn 



D 
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where w(x,y) = E[t\x(0) = x,y(0) = y] satisfies A( Xj y\W = —1/D. The 
Laplacian transforms according to 

&(x,v)W = (v>l + u 2 y )A {U)V) w, 

by the Cauchy-Riemann equations and the Jacobian of the transfor- 
mation is J = u 2 + u 2 . Therefore, 

j A( x ^wdxdy= / A^^w dudv. 

This means that the compatibility condition of Section [2] remains 
unchanged and gives the area of the original domain. 



5 Domains with cusps 

Here we find the leading order term of the MFPT for small holes located near 
a cusp of the boundary. A cusp is a singular point of the boundary. As a = 
at the cusp, one expects to find a different asymptotic expansion than (|4.1(Jj) . 
As an example, consider the Brownian motion inside the domain bounded 
between the circles (x — 1/2) 2 + y 2 = 1/4 and (x — 1/4) 2 + y 2 = 1/16 (see 
Figure EJ). The conformal mapping z i— ► exp{iri(l/z — 1)} maps this domain 
onto the upper half plane. Therefore, the MFPT is to leading order 

ET=§(l +0m ). (5.1) 

This result can also be obtained by mapping the cusped domain 
to the unit circle. The absorbing boundary is then transformed to 
an exponentially small arc of length exp{— rr/e} + 0(exp{— 27r/e}), and 
equation (J5.1|) is recovered. 

If the ratio between the two radii is d < 1, then the conformal map 
that maps the domain between the two circles to the upper half plane is 

{7TZ I 
— [\jz — 1) > (for d = 1/2 we arrive at the previous example), 
d' l -l J 

so the MFPT is to leading order 

The MFPT tends algebraically fast to infinity, much faster than the O (log - 

V £ 

behavior near smooth or corner boundaries. The MFPT for a cusp is much 
larger because it is more difficult for the Brownian motion to enter the cusp 
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than to enter a corner. The MFPT (15. 2 p can be written in terms of d 
instead of the area. Substituting |Q| = 7ri? 2 (l — d 2 ), we find 

^^§±4(I + 0(1) ), ,5.3, 

where R is the radius of the outer circle. Note that although the 
area of is a monotonically decreasing function of d, the MFPT is 
a monotonically increasing function of d and tends to a finite limit 

as d — > 1. 

Similarly, one can consider different types of cusps and find that the leading 
order term for the MFPT is proportional to l/e x , where A is a parameter that 
describes the order of the cusp, and can be obtained by the same technique of 
conformal mapping. 



6 Diffusion on a 3-sphere 
6.1 Small absorbing cap 

Consider a Brownian motion on the surface of a 3-sphere of radius R 
described by the spherical coordinates (9, 0) 

x = Rsin9 cos0, y = R sin 9 sin 0, z = Rcos8. 

The particle is absorbed when it reaches a small spherical cap. We center the 
cap at the north pole, 9 = 0. Furthermore, the FPT to hit the spherical cap 
is independent of the initial angle 0, due to rotational symmetry. Let v{9) be 
the MFPT to hit the spherical cap. Then v satisfies 

A M v = -l, (6.1) 

where Am is the Laplace-Beltrami operator ^H] of the 3-sphere. This Laplace- 
Beltrami operator Am replaces the regular plane Laplacian, because the dif- 
fusion occurs on a manifold [HH and reference therein]. For a function v in- 
dependent of the angle the Laplace-Beltrami operator is (see Appendix IX]) 

A M v = R- 2 (v" + cot 9 v') . (6.2) 
The MFPT also satisfies the boundary conditions 

v'(n) = 0, v{5) = 0, (6.3) 

where 5 is the opening angle of the spherical cap. The solution of the boundary 
value problem ()6.2|) . (|6.3j) is given by 

vm = 2RHo g ^m y (6.4) 
sm(<)/2) 
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Not surprisingly, the maximum of the MFPT is attained at the point 6 = it 
with the value 



vmax = v{ir) = -2R 2 log sin 5 - = 2R 2 (bg ^ + log 2 + 0(S 2 )^j . (6.5) 
The MFPT, averaged with respect to a uniform initial distribution, is 

1 r 

Et = j / v(6)sm6d0 

2 cos 2 - Js 
2 

2 f \og S m(5/2) , 1\ 
V cos 2 (<5/2) 2 J 

= 2R 2 ^log^ + log2-i + 0(5 2 log5)^ . (6.6) 
Both the average MFPT and the maximum MFPT are 

\n\ n / i 



\<j 



log - + 0(1) , (6.7) 



2tt V 5 

where \Q\ g = AttR 2 is the area of the 3-sphere. This asymptotic expansion 
is the same as for the planar problem of an absorbing circle in a disk. The 
result is two times smaller than the result (|1.1|) that holds when the absorbing 
boundary is a small window of a reflecting boundary. The factor two difference 
is explained by the aspect angle that the particle "sees". The two problems 
also differ in that the "narrow escape" solution is almost constant and has a 
boundary layer near the window, with singular fluxes near the edges, whereas 
in the problem of puncture hole inside a domain the flux is regular and there 
is no boundary layer (the solution is simply obtained by solving the ODE). 

6.2 Mapping of the Riemann sphere 

We present a different approach for calculating the MFPT for the Brownian 
particle diffusing on a sphere. We may assume that the radius of the sphere is 
1/2, and use the stereographic projection that maps the sphere into the plane 
[T7j . The point Q = (£,rj, () on the sphere (often called the Riemann sphere) 

£ 2 + r? 2 + (C-l/2) 2 =(l/2) 2 
is projected to a plane point P = (x, y, 0) by the mapping 
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and conversely 



t = TTTi> 7 l = TT-2' C = TT^lr (6-9) 
1 + r 2 1 + r l 1 + r l 

The stereographic projection is conformal and therefore transforms harmonic 
functions on the sphere harmonic functions in the plane, and vice versa. How- 
ever, the stereographic projection is not an isometry. The Laplace-Beltrami 
operator Am on the sphere is mapped onto the operator (1 + r 2 ) 2 A in the 
plane (A is the Cartesian Laplacian). The decapitated sphere is mapped onto 
the interior of a circle of radius 

r 

r* = cot-. (6.10) 

Therefore, the problem for the MFPT on the sphere is transformed into the 
planar Poisson radial problem 

^V=- 1 for r<r s , (6.11) 

(1 + r A y 

subject to the absorbing boundary condition 

V(r = r 5 ) = 0, (6.12) 

where 

V(r) = v{6). 

The solution of this problem is 

Transforming back to the coordinates on the sphere, we get 

, , 1, sin(0/2) 
2 sin(o/2J 

As the actual radius of the sphere is R rather than 1/2, multiplying eq. ()6.14|) 
by (2i?) 2 , we find that (l6~T4l is exactly 

6.3 Small cap with an absorbing arc 

Consider again a Brownian particle diffusing on a decapitated 3-sphere of ra- 
dius 1/2. The boundary of the spherical cap is reflecting but for a small window 
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that is absorbing (see Fig. [H}. We calculate the mean time to absorption. Us- 
ing the stereographic projection of the preceding subsection, we obtain the 
mixed boundary value problem 

Av = -— — for r <r s , < 4> < 2n, 



v(r, <p) 
dv(r, (p) 



dr 

The function 



r=r s 



r=r s 



(1 + r 2 
0, for 

0, for 



7r < e, 



7T > £. 



(6.15) 



1 , ( 1 + rj 



w (r) = - log 

I y i _|_ r ^ 



is the solution of the all absorbing boundary problem eq. (|6.13|) . so the function 
u = v — w satisfies the mixed boundary value problem 



Am 



0, r < r 5) for < < 2tc, 



u(r, <p) 
du(r, <fi) 



dr 



r=r s 



r=r s 



0, for \(j> — tt\ < e, 
rs 



2(1 + rl) 



2\> 



for 



7r > e. 



(6.16) 



Scaling f = r /rs, we find this mixed boundary value problem to be that of a 
planar disk [S], with the only difference that the constant 1/2 is now replaced 



by 



2(1 + r? 



-. Therefore, the solution is given by 



a 



2rf 



log - + 0(e) 



(6.17) 



Transforming back to the spherical coordinate system, the MFPT is 
v{d,<f>) 



1 sin^/2 

2 ° g ^J/2 



(6.18) 



5 



cos 



n=l 



cot (0/2) 
cot (5/2) 



cos nc 



The MFPT, averaged over uniformly distributed initial conditions on the de- 
capitated sphere, is 



Et = — 
2 



1 /logsin(5/2) 1 



cos 2 (5/2) 



• - I +cos 2 i 



log - + 0(e) 



(6.19) 
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Scaling the radius R of the sphere into (|6.19|) . we find that for small e and S 
the averaged MFPT is 



Et = 2R 2 



logi + 21og- + 31og2 

d e 



+ 0(e,6 2 log 5, S 2 \oge) 



(6.20) 



There are two different contributions to the MFPT. The ratio e between the 
absorbing arc and the entire boundary brings in a logarithmic contribution to 
the MFPT, which is to leading order 



\n\ 



7T 



log- 



However, the central angle 5 gives an additional logarithmic contribution, of 
the form 



2tt 



logi. 



The factor 2 difference in the asymptotic expansions is the same as encountered 
in the planar annulus problem. 

The MFPT for a particle initiated at the south pole 9 = it is 



v(tt) 



S S 
-2R 2 log sin - - AR 2 cos 2 - 
6 2 2 



log - + 0(e) 



2R 2 



log T + 2 log - + 3 log 2 + O (e, 5 2 log 5, 5 2 log e 

o e 



. (6.21) 

for which the MFPT is maximal. The sta- 



We also find the location (6, 

dv 

tionarity condition — — = implies that 6 = 0, as expected (the opposite 

^-direction to the center of the window). The infinite sum in equation (j6.18|) 
is 0(1). Therefore, for 5 <ti 1, the MFPT is maximal near the south pole 
9 = tt. However, for 5 = 0(1), the location of the maximal MFPT is more 
complex. 

Finally, we remark that the stereographic projection also leads to the de- 
termination of the MFPT for diffusion on a 3-sphere with a small hole as 
discussed above, and an all reflecting spherical cap at the south pole. In this 
case, the image for the stereographic projection is the annulus, a problem 
solved in Sectional 



A Laplace Beltrami operator on 3-sphere 

The Laplace Beltrami operator on a manifold is given by 

Am/ = -=L= £ J- ( 9 «VteG°£) , 
VdetG d& \ d^J 
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where 

ti = Wi' ^' = ^'*^' G = ^' tf 3 =^- ( A - 2 ) 
In spherical coordinates we have 

g ee = R 2 , 944 = R 2 sin 2 9, g <t> = g<po = 0. (A. 3) 

Therefore, for a function w — w(0, <fi) 

AM/ .^(0 + cot9 | + _i_p). (A . 4) 
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Figure 1: An annulus R\ < r < R 2 . The particle is absorbed at an arc of 
length 2eRi (dashed line) at the inner circle. The solid lines indicate reflecting 
boundaries. 
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Figure 2: Rectangle of sizes a and b with a small absorbing segment of size 
at the corner. 
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e 

•4 ► 




Figure 3: A small opening near a corner of angle a. 
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Figure 4: The point (0, 0) is a cusp point of the dotted domain bounded 
between the two circles. The small absorbing arc of length e is located at the 
cusp point. 



Figure 5: A sphere of radius R without a spherical cap at the north pole of 
central angle 5. The particle can exit through an arc seen at angle 2e. 
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